Draft version February 2, 2008 

Preprint typeset using LAT^X style emulatcapj v. 04/03/99 



A NOTE ON GRAVITATIONAL BROWNIAN MOTION 
David Merritt 

Department of Physics, Rochester Institute of Technology, Rochester, NY 14623 drmsps@ad.rit.edu 

Draft version February 2, 2008 

ABSTRACT 

Chandrasekhar's theory of stellar encounters predicts a dependence of the Brownian motion of a mas- 
r^j- ■ sive particle on the velocity distribution of the perturbing stars. One consequence is that the expectation 

value of the massive object's kinetic energy can be different from that of the perturbers. This effect is 
shown to be modest however, and substantially smaller than claimed in a recent study based on an 
£SJ ■ approximate treatment of the encounter equations. 

Subject headings: black hole physics — gravitation — gravitational waves — galaxies: nuclei 
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1. INTRODUCTION 



A massive object at the center of a stellar system undergoes a random walk in momentum space as its motion is perturbed 
by gravitational encounters with nearby stars. After the decay of transients, the steady-state velocity distribution f(v) 
of the massive object is a balance between the accelerating forces of random perturbations and the decelerating force of 
[ dynamical friction. Chandrasekhar's (1942,1943) theory of gravitational encounters provides expressions for both forces 
■ as functions of the mass and velocity of the test particle, and the masses and velocity distribution of the field stars. In the 
j— ( ' case that the latter have a Maxwellian velocity distribution, the steady-state f(v) describing the motion of the massive 
I/"") , object is also a Maxwellian, with rms velocity given by the equipartition relation 
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, where M and v refer to the massive object and m and u to the perturbers 

One difference between Brownian motion in gravitating and non-gravitating systems is that the velocity distribution 

//(it) of the perturbers in a galaxy or star cluster can be substantially non-Maxwellian. The steady-state velocity 

distribution of the massive object is still predicted to be Maxwellian (Merritt 2001), but there is no reason for v rms to 

O ' satisfy equation QJ. In this note, the steady-state f(v) is derived for a massive particle that is subject to gravitational 

perturbations from field stars with an arbitrary distribution of velocities. Departures from equipartition are found to be 

cZ) ■ small, even when the background is strongly non-Maxwellian. The stronger dependence of v rms on ff(u) claimed in a 

' recent study (Chatterjee, Hernquist & Loeb 2002, hereafter CHL) is shown to result from approximations made in the 

, application of Chandrasekhar's theory. 
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^ ■ 2. STEADY-STATE VELOCITY DISTRIBUTIONS 

The steady-state velocity distribution f(v) of a massive object ("black hole") that interacts via gravitational scattering 
with stars is given by the time-independent solution of the Fokker-Planck equation. Following the treatment in Merritt 
(2001) (hereafter Paper I), we can write 
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where (Avy), (Aujj) and (Aw^) are the standard diffusion coefficients describing the motion of the black hole. In a steady 

state, we expect that v will be of order yJm/M <C 1 times the typical stellar velocity. Expanding the diffusion coefficients 
about v = 0, 

(Au||) = -Av + Bv 3 . . . , (3a) 

(A«jf) = C + Dv 2 . . . , (3b) 

(Avl) =2(C + Fv 2 )..., (3c) 

substituting into equation ©, and retaining the lowest order terms in v gives 

0^(A + 2D-F)f + ^- (4a) 
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since A is of order M/m times D and F. This has solution 
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i.e. the black hole's velocity follows a Maxwellian distribution with ID velocity dispersion a. 

If the background stellar distribution is assumed to be infinite and homogeneous, the diffusion coefficients are (Paper 



I) 
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(Aujj) = —ir 2 G 2 m 2 nv / du \~) ff(u)H 2 (v, u,p max ), 
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where ff(u) is the velocity distribution of the field stars, normalized to unit number; n is the field star density; p max is 
the maximum impact parameter at which incoming stars are assumed to perturb the black hole; and 
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Equation (|7b|) includes the "non-dominant" terms which are of order unity when p ma x is large. The integration variable 
V is the relative velocity of star and black hole when the star is at infinity. 

In the standard approximation (e.g. Rosenbluth, MacDonald & Judd 1957; Spitzer 1987), the non-dominant terms are 
neglected, and the logarithmic terms are taken outside the integrals, since they are slowly varying with respect to V. One 
writes 
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a constant, and 
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with leading terms 
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A A = —Tr 2 G 2 MmnlnAf f (0), 
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The predicted ID velocity dispersion of the black hole is then (CHL) 
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However when the velocity of the massive particle is very low, the logarithms in equations (JJaJ-b) will be close to zero 
for all stars with u < v, i.e. all stars that enter into the dynamical friction integral l|9a|l . It is unclear in this case whether 
removing the logarithm from the integrands is a reasonable approximation (Chandrasekhar 1943; White 1949). 

Returning to the exact expressions (|6al7b|l and expanding to lowest order in v, we find after some algebra: 
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The non-dominant terms may be shown not to contribute to this order. Note that both diffusion coefficients now depend 
on the field-star velocity distribution at all values of u, and in fact the field stars that produce the dynamical friction force 
can be shown to have roughly the same velocity distribution as ffiu) (Paper I). In the standard treatment with "logA's 
removed," the dynamical friction force comes entirely from stars with u < v (as in eauation l9af> . a physically unreasonable 
result. 

The black hole's velocity dispersion is now 

2 _ C _ 1 m io°° duu ff( u ) lo s (i + 2 m^rr) 
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If ff(u) is differentiable at low and high u, the denominator can be integrated by parts yielding 
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This reduces to the approximate expression, equation H 1 1 p . if the logarithmic terms are "taken out of the integrals." 
In the case of a Maxwellian field star velocity distribution, //(it) = fh[{u) — foe~ u ^ 2(Jf , the exact coefficients are 

-Am = — t, % G(Rm), Cm = — 7, G(Rm), (15a) 
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and 
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independent of p ma x (Paper I) . Remarkably, this is also the value of a 2 given by the approximate expression (|11|) when 
//(«) = Im{u). 

For all other - non-Maxwellian - //'s, however, equations l|13|l and (|11() yield different results for cr 2 , and each is 
different from the Maxwellian value. Furthermore a 2 as predicted by the exact expression (|13fl is in general a function of 
the maximum impact parameter, since C and A depend differently on p m ax- 

Without loss of generality, we can express the dependence of a 2 on p max y i a a dimensionless variable i?, where R = 
PmaxO' 2 /GM and a is a ID velocity dispersion that characterizes the field star velocity distribution. Equating G(R) in 
equation l|15b(l with In A = ln(p m ax/Pmin) reveals that p m in ~ GM / 1 \[2a y (Paper I), hence R w Pmax/Pmin'i note that 
Pmin is also approximately equal to the "radius of influence" of the black hole. When p max is large, i.e. R 3> 1, the 
logarithmic terms in equation (|14|) are nearly constant and a 2 reduces to the value (|11|) derived by "taking out In A." But 
for a supermassive black hole in a galactic nucleus, p m ax is of order the core radius r c of the stellar distribution (Paper 
I), and observed core radii are of order a few times GM/a 2 at most (e.g. Poon & Merritt 2001); hence p ma x > Pmin and 
R is of order unity. Thus we expect the exact and approximate expressions for a 2 to be substantially different in many 
astrophysically interesting situations. 

3. EXAMPLES 

The equations derived above can be applied to various field-star distributions and the results for a compared with the 
equipartition value. Following CHL, we quantify the departure of a from equipartition via the parameter 77, defined as 

- Mcj2 (^7^ 
to a+ 

with cr* the stellar velocity dispersion at the location of the black hole; 77 is unity in the case of equipartition. Both a and 
cr* are computed under the assumption that the black hole sits at the center of the stellar system. We emphasize that 
the effect of the black hole on the stellar background is ignored. 
King's (1966) family of models 

me) = { (^r 3/2 - 1) ■ if ^°; as) 

[ 0, otherwise 

have nearly Maxwellian central velocity distributions. Here E is the binding energy, E = — u 2 /2 + ip(r), and ip(r) = — $(r) 
with $(r) the gravitational potential due to the stars, equal to zero at the model's edge. The parameter W = ^{Q)/a 2 K is 
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a dimensionless measure of the central potential. For large W, the King model parameter <jk is nearly equal to the true 
central velocity dispersion o+, or 

j2 _2 CdE{W-Ef 2 f K {E) 
* 3 dEVW=Ef K (E) 

A natural definition for R is 

2 

Rk = ~gnT- (20) 

King models having W ^> 10 fit the observed brightness profiles of some elliptical galaxies moderately well. Setting M to 
~ 0.1% of the total stellar mass and equating p max to the core radius gives Rk of order unity when W ^ 10. Figure 1 
shows that r\ is extremely close to one when R and W have these values. King models with smaller W do not represent 
galaxies well - their cores are too large - but are good fits to some globular clusters. The appropriate value of Rk for 
globular clusters is unclear however since it is not known whether such systems harbor massive central objects or what 
their likely masses are. If we suppose that the black hole has a mass equal to a fraction F of the stellar mass in the core, 
then Rk ~ F 1 " 1 . An upper limit on R comes from assuming that the black hole has a minimum mass ~ IOMq and that 
the core mass is ~ 5 x 10 4 M©, its value in a massive, large-cored globular cluster like Li Centauri (e.g. Merritt, Meylan 
& Mayor 1997). Then Fk2x 10~ 4 and Rk ~ 5000. In such an extreme case, Figure 1 shows that rj can approach its 
asymptotic limit of 1.75 (CHL), but only in highly unphysical models with W < 2. 

Next we consider the polytropic models, some of which have strongly non-Maxwellian velocity distributions. The 
Plummer (1911) polytrope has // = /5(E), where 

fn(E) = \f° En - l > *jf>?5 (21) 

0, otherwise 

and the gravitational potential is again defined to be zero at the model's edge, which in this case lies at infinity. The 
density profile is 

3M gal a 2 1 

P\ r ) = 7 fT77 ( 22 ) 
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and the central velocity dispersion is <r 2 (0) = GM ga i/6a. We define 
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which gives 



Rp = GM (23) 



= 6 ^(2-x)^log(l + 36i^ 2 ) 

V 7 r 2 (fe(2-x) 5/2 log(l + 36i? 2 a; 2 )' 

As R — > co, rj — > 4/3 (CHL). Figure 2 shows i] as a function of R. Values of R appropriate to supermassive black holes 
in galactic nuclei again give rj close to one and the asymptotic value is approached only slowly as R is increased. 

The third example considered by CHL was the n = 1 polytrope. The structure of this model makes it a poor represen- 
tation of any stellar system, since the density profile: 

P( r ) = . T\ r ~ 1 sin(r/o), r < ira (25) 

is "all core." Furthermore the distribution function is singular at the model's edge. For this reason we must use the first 
of the two expressions derived above for a 2 , equation (|13|) . since the second assumes differentiability of //. We define R 
in terms of the central velocity dispersion as in equation l|23l) . where cr 2 (0) = GM ga i/2ira. The result for r\ is 

_ = _J_ ^d X ^- x y 1/2 \o g (i + i^R 2 ) 

Figure 3 plots this function. The asymptotic value r\ = 4 (CHL) is again reached only very gradually as R increases. 

4. DISCUSSION 

The foregoing should not be taken to imply that a massive object at the center of a stellar system should always be in or 
near a state of equipartition with respect to its less-massive perturbers. Indeed it is not even clear what Chandrasekhar's 
theory would predict in many physically interesting situations. Chandrasekhar's theory is local in the sense that the 
density and velocity distribution of the field stars are assumed to everywhere be the same as they are at the location of 
the test particle. The effect of a massive particle on the structure of the background stellar system is also ignored. Both 
of these assumptions make it difficult to apply the theory to the case of Brownian motion of a supermassive black hole 
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at the center of a galactic nucleus. The density of stars falls steeply away from the black hole in most galaxies; indeed a 
density falloff as steep as p ~ r -1 / 2 is required for self-consistency if the stellar velocities are isotropic near the black hole 
(as assumed here). Furthermore the presence of the black hole implies a steeply-rising stellar velocity dispersion within 
its sphere of influence, r <; GM/a%, and the black hole's Brownian motion would be influenced by these higher velocities. 

The likely importance of such effects is suggested by a recent numerical study of Brownian motion in a galaxy with a 
steep, p <~ r~ 3//2 density cusp (Dorband, Hemsendorf & Merritt 2003). The stellar velocity dispersion in this model has a 
central value of zero in the absence of the black hole, yet the numerical integrations reveal a significant Brownian motion 
of the "black hole." Evidently the massive particles in these simulations are responding to stars far from the center, or to 
stars whose motions have themselves been influenced by the presence of the black hole, or both. 

This work was supported by grants AST-0071099 and AST-0206031 from the NSF and by grant MAG5-9046 from 
NASA. 
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Fig. 2. — Values of r] f or a massive object at the center of a Plummer-model galaxy, as a function of the dimensionless maximum impact 
parameter R of equation 1231 . Upper dashed line shows asymptotic value for large R. 




Fig. 3.— 



Like Figure 2, for the n = 1 polytropc. 



